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This paper proposes enhanced fruit fly optimization algorithm (EFF) 
and status of material algorithm (SMA) to solve the optimal reactive power 
problem. Fruit fly optimization algorithm is based on the food finding 
behavior of the fruit fly. There are two steps in food finding procedure of 
fruit fly: At first it smells the food source by means of osphresis organ and it 
flies in that direction; afterwards, when it gets closer to the food site, through 
its sensitive vision it will find the food. At the beginning of the run by 
diminishing the inertia weight from a large value to a small value, will lead 
to enhance the global search capability and more local search ability will be 
in process the end of the run of the EFF algorithm. Then SMA is projected to 
solve the problem. Three state of material are solid, liquid, and gas. For 
evolution procedure direction vector operator assign a direction to every 
molecule consecutively to guide the particle progression. Collision operator 
imitates the collisions factor in which molecules are interacting to each other. 
Proposed enhanced EFF, SMA has been tested in standard IEEE 30 bus test 
system and simulation results show the projected algorithms reduced the real 
power loss considerably. 
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1. INTRODUCTION 

Reactive power problem plays an important role in secure and economic operations of power 
system. Numerous types of methods [1-6] have been utilized to solve the optimal reactive power problem. 
However many scientific difficulties are found while solving problem due to an assortment of constraints. 
Evolutionary techniques [7-15] are applied to solve the reactive power problem. This paper proposes 
enhanced fruit fly optimization algorithm, status of material algorithm to solve the optimal reactive power 
problem. Fruit fly optimization algorithm is based on the food finding behavior of the fruit fly. When 
compare to other species fruit fly is better to in vision and osphresis. In the projected enhanced fruit fly 
optimization algorithm (EFF), linear generation mechanism of candidate solution is incorporated with fruit 
fly algorithm. In order to maintain the balance between global and local search an inertia weight is applied 
in the procedure. For global search large inertia weight is utilized and small inertia weight applied for local 
search. At the beginning of the run by diminishing the inertia weight from a large value to a small value, 
will lead to enhance the global search capability and more local search ability will be in process the end 
of the run. Then status of material algorithm (SMA) is utilized to solve the optimal reactive power problem. 
Three state of material are solid, liquid, and gas. In gas there is no precise shape, but fit into the entire 
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container in which it is restricted. Development is tested by the molecules represent the utmost allowable 
displacement pi amongst the particles. Intermolecular forces restricted in liquid state, than in gas state. 
There will be sufficient energy for the molecules to shift reasonably to each other but a mobile structure will 
reign. Throughout the container the shape of the liquid can be accessed. Molecules denote a particle 
movement p2 within the liquid state. Molecules are crammed collectively with forces amongst the particles 
are being tough enough cannot move liberally only can vibrate in the solid state. Solid state has a steady, 
unambiguous shape and possesses specific volume. Solids can only modify their shape by applying force, 
and particles are capable to vibrate taking into account of a minimum p3 distance. Collision operator imitates 
the collisions factor in which molecules are interacting to each other. In arbitrary behaviour of molecules, 
capricious positions are formed subsequent to a probabilistic condition which considers capricious locations 
within a realistic exploration space. Proposed enhanced fruit fly optimization algorithm (EFF), status 
of material algorithm (SMA) has been tested in standard IEEE 30 bus test system and simulation results show 
the projected algorithms reduced the real power loss considerably. 


2. PROBLEM FORMULATION 

Objective of the problem is to reduce the true power loss: 


F = P L = IkeNbr gk (^ 2 + Vj 2 - ZViVjCOSOij) 

(1) 

Voltage deviation given as follows: 

F = P L + cl> v x Voltage Deviation 

(2) 

Voltage deviation given by: 

Voltage deviation = Sjl^lVj — 1| 

(3) 

Constraint (equality): 

Pg — Pd + Pl 

(4) 

Constraints (inequality): 

pmin p pmax 

r gslack — r gslack — r gslack 

(5) 

<Q gi <Q™ ax ,i£N g 

(6) 

Vtm" < v ; < Vi max , i £ N 

(7) 

T min < T . < T .max ; j £ Nt 

(8) 

Q™ in < Qc < Qc ax , i £ N c 

(9) 


3. ENHANCED FRUIT FLY OPTIMIZATION ALGORITHM 

Fruit fly optimization algorithm is based on the food finding behavior of the fruit fly. When 
compare to other species fruit fly is better to in vision and osphresis [16]. There are two steps in food finding 
procedure of fruit fly: at first it smells the food source by means of osphresis organ and it flies in that 
direction; afterwards, when it gets closer to the food site, through its sensitive vision it will find the food. 
Modeled algorithm of the fruit fly optimization has been given below: 

a. Initialization of parameters 

b. Candidate solutions are engendered 

c. Location of fruit fly swarm at preliminary level 

x_axis = random ( LR ) ; y_axis — random ( LR ) (10) 

d. Confer capricious direction and distance for finding the food by an individual fruit fly 
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Xj = x_axis + random ( FR ) ; y t = y_axis + random (FR) 

e. From the origin compute the distance of food location 
Distance t = jxf + yf 


f. Value of the smell concentration judgment is found by. 


S t = —- 

Distance, 


Si = 


1 J(_x_axis +random ( FR)) 2 +(y_axis +random ( FR ) ) 2 


(ID 


( 12 ) 


(13) 

(14) 


g. Individual fruit fly smell concentration is calculated by, 

Smelli = objectivefunction(Si) (15) 

h. Among the fruit fly swarm discover out the fruit fly with maximum smell concentration, 

[smell best, best index] — maximum(smell ) (16) 


i. By using vision towards that location fruit fly swarm flies, 

smell best = best smell (17) 

x_axis = x(best index ) (18) 

y_axis — y(best index ) (19) 

j. Replicate the implementation of steps b-i. When the smell concentration is not superior to 
the preceding iterative smell concentration any longer or when the iterative number reaches 
the maximum number of iterations, the flow will stop. 

In the projected enhanced fruit fly optimization algorithm (EFF), linear generation mechanism 
of candidate solution is incorporated with fruit fly algorithm. In order to maintain the balance between global 
and local search an inertia weight is applied in the procedure. For global search large inertia weight is utilized 
and small inertia weight applied for local search. At the beginning of the run by diminishing the inertia 
weight from a large value to a small value, will lead to enhance the global search capability and more local 
search ability will be in process the end of the run. 

a. Initialization of parameters 

b. Candidate solutions are engendered through linear generation mechanism 

c. Location of fruit fly swarm at preliminary level 

x_axis' — n * random (20) 


d. Confer capricious direction and distance for finding the food by an individual fruit fly 

x'i = x_axis' + a) * random \a> = a> 0 * a a eneratwn (21) 

e. Value of the smell concentration judgment is found by 

S[ = x[ — xjaxis' + a) * random (22) 

k. Individual fruit fly smell concentration is calculated by. 

Smell • = objectivefunction{S[) (23) 

l. Among the fruit fly swarm discover out the fruit fly with maximum smell concentration, 

[smell best', best index'] = maximum(smell ') (24) 
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m. Maintain the maximum concentration value and ^'coordinate. subsequently, by using vision the fruit 
fly swarm flies in the direction of that location: 

smell best’ — best smell' (25) 

x_axis' = x(best index)' (26) 

n. Replicate the implementation of Steps b-i. When the smell concentration is not superior to 
the preceding iterative smell concentration any longer or when the iterative number reaches 
the maximum number of iterations, the flow will stop. 


4. STATUS OF MATERIAL ALGORITHM 

Conventionally, three state of material are solid, liquid, and gas. In gas there is no specific shape, 
but fit into the whole container in which it is restricted. Progression is tested by the molecules symbolize the 
utmost allowable displacement pi amongst the particles. Intermolecular forces are extra restricted in liquid 
state, than in gas state. There will be enough energy for the molecules to shift comparatively to each other but 
a mobile structure will prevail. Through the container the shape of the liquid can be accessed. Molecules 
symbolize a particle movement p2 inside the liquid state. Molecules are crammed collectively with forces 
amongst the particles are being tough enough cannot move liberally only can vibrate in the solid state. Solid 
state has a steady, unambiguous shape and possesses specific volume. Solids can only modify their shape by 
applying force, and particles are capable to vibrate taking into account of a minimum p3 distance [17], 

In the proposed method status of material algorithm (SMA), for evolution procedure direction vector 
operator assign a direction to every molecule consecutively to guide the particle progression. Collision 
operator imitates the collisions factor in which molecules are interacting to each other. In arbitrary behaviour 
of molecules, capricious positions are formed subsequent to a probabilistic condition which considers 
capricious locations within a realistic exploration space. Primarily within the range of [-1, 1], all the direction 
vectors (c = are arbitrarily chosen. Molecules experience numerous attraction forces 

and new-fangled direction vector is calculated by: 


9i +1 — Si (l — 0.5 + a t 

(27) 

a l = (Q best -Q i )/\\Q best -Q l \\, 

(28) 

Velocity V t of every molecule is calculated by. 


~ Si' Vinit 

(29) 

Initial velocity v initial magnitude computed by. 


v n (u hi 9 h ulow\ 

_ bj ) 

v initial n P 

(30) 

New-fangled location for every molecule is modernized by, 


Qij 1 = Qtj + vtj ■ rand( 0,1) ■ p ■ ( b bigh - b) ow ) 

(31) 

When two molecules which is shorter than a determined proximity value then Collisions will occur. 
Consequently, if ||<2; - f) ( J| < r, a collision between molecules i and q is assumed; or else, there won’t be 
collision, in view of i,q E N P ] -t q. When collision occurs, every particle direction vector is tailored by 

swapping their particular direction vectors as follows: 

Si Sq ttnd Sq Si 

Collision radius computed by. 

(32) 

n 

(33) 
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Within the range [0, 1] an uniform arbitrary number r,„is generated and modelled by, 

Qk+i | bj ow + rancL( 0,1) ■ ( b* l0h — bj ow ) with probability H 
( Qwith porbability 1 — H 

Most excellent found individual from the current k population Q htsl ’ k is compare to the most 
excellent individual Q bt "’ l ' k ' 1 of the previous generation. If Q besf - k is enhanced than Q bcst ' k ‘' by fitness value, 
then best Q is modernized with Q best-k , or else Q best remains without any change. Consequently, Q best stock up 
the most excellent individual found up to now. Modus operandi map the existing population Q k to the new-fangled 
population Q ,l+I . Algorithm obtain the existing population Q A and the configuration parameters p, (3, a and H, 
yield the new-fangled population Q* +l . 

- Common procedure: 

Step a: Best element of the population Q best £ {Q } has been found 
Step b: venial, r are calculated 

Step c: By means of the Direction vector operator find the new-fangled molecules 
Step d: By utilizing the Collision operator resolve the collisions 

Step e: By utilizing the arbitrary positions operator find the new-fangled arbitrary positions 

In the optimization procedure, 50% of iterations for the gas state (exploration), 40% of iterations 
for the liquid state (exploration-exploitation) and 10% of iterations for the solid state (exploitation). 
Algorithm begins by initializing a set Q of Nq molecule (q = {&, Q 2 ,.., each molecule position Q; is 
an-dimensional vector containing the parameter values to be optimized. Values are arbitrarily 
and unvaryingly dispersed between the pre-specified lower preliminary parameter bound bj ow and the upper 
initial parameter bound b^ L ‘ qh as follows, 

Qfj = bj ow + rand (0,1) ■ (bf ah - b) ow ) (35) 

- Gas status 

Step a: Parameters p £ [0.80, 1], /?=0.80, a=0.80, and H= 0.90 are fixed 
Step b: Common procedure will be applied 

Step c: When 50% of the entire iteration number is ended (1 <£<0.5 < generation), afterwards the procedure 
continues to the liquid state process; or else go back to step b. 

- Liquid status 

Step d: Parameters p £[0.30, 0.60], /?=0.40, a=0.20, and H= 0.20 are fixed 
Step e: Common procedure will be applied 

Step f: When 90% (50% from the gas state and 40% from the liquid state) of the entire iteration number is 
concluded (0.5 -generation<k< 0.9 ■ generation ), afterwards the procedure continue to the solid state process; 
if not go back to step e. 

- Solid status 

Step g: Parameters p £[0.0, 0.1] and /?=(). 1, a=0, and H= 0 are fixed 
Step h: Common procedure will be applied 

Step i: If the 100% of the entire iteration number is concluded (0.9 -gen<k<gen), then the procedure is 
completed; if not go back to step h. 


5. SIMULATION RESULTS 

Proposed enhanced fruit fly optimization algorithm (EFF) and status of material algorithm (SMA) 
has been tested in standard IEEE 30 Bus system [18]. Table 1 shows the constraints of control variables. 
Table 2 shows the limits of reactive power generators and comparison results are presented in Table 3. 
Figure 1 gives the comparison of real power loss and Figure 2 gives the reduction of real power loss (%) with 
reference to base case value. 


Table 1. Constraint values of the control variables 


System type 

List of variables 

Minimum value (PU) 

Maximum value (PU) 

IF.F.F. 30 Bus 

Generator voltage 

0.95 

i.i 


Transformer tap 

0.9 

i.i 
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VAR source 

0 

0.20 

Table 2. Constrains values of the reactive power generators 

System type 

List of variables 

Value of Q minimum (PU) 

Value of Q maximum (PU) 

TF.F.F, 30 Bus 

i 

0 

10 


2 

-40 

50 


5 

-40 

40 


8 

-10 

40 


11 

-6 

24 


13 

-6 

24 


Table 3. Simulation results of IEEE-30 system 


List of control variables 

Base case value 

MPSO [19] 

PSO [19] 

EP [19] 

SARGA [19] 

EFF 

SMA 

VG-l 

1.060 

1.101 

1.100 

NR* 

NR* 

1.014 

1.021 

VG-2 

1.045 

1.086 

1.072 

1.097 

1.094 

1.030 

1.029 

VG-5 

1.010 

1.047 

1.038 

1.049 

1.053 

1.011 

1.010 

VG-8 

1.010 

1.057 

1.048 

1.033 

1.059 

1.019 

1.018 

VG-l 2 

1.082 

1.048 

1.058 

1.092 

1.099 

1.020 

1.016 

VG-l 3 

1.071 

1.068 

1.080 

1.091 

1.099 

1.024 

1.020 

Tap-11 

0.978 

0.983 

0.987 

1.01 

0.99 

0.942 

0.937 

Tap-12 

0.969 

1.023 

1.015 

1.03 

1.03 

0.911 

0.921 

Tap-15 

0.932 

1.020 

1.020 

1.07 

0.98 

0.909 

0.914 

Tap-36 

0.968 

0.988 

1.012 

0.99 

0.96 

0.917 

0.932 

QC-10 

0.19 

0.077 

0.077 

0.19 

0.19 

0.090 

0.090 

QC-24 

0.043 

0.119 

0.128 

0.04 

0.04 

0.105 

0.104 

PG (MW) 

300.9 

299.54 

299.54 

NR* 

NR* 

298.37 

298.41 

QG (Mvar) 

133.9 

130.83 

130.94 

NR* 

NR* 

130.18 

130.26 

Reduction in PLoss (%) 

0 

8.4 

7.4 

6.6 

8.3 

18.23 

20.38 

Total PLoss (Mw) 

17.55 

16.07 

16.25 

16.38 

16.09 

14.349 

13.972 


Note: NR*-Not reported 




Figure 1. Comparison of real power loss Figure 2. Reduction of real power loss (%) with 

reference to base case value 


6. CONCLUSION 

Proposed enhanced fruit fly optimization algorithm (EFF), status of material algorithm (SMA) 
successfully solved the optimal reactive power problem. In the projected Enhanced Fruit fly optimization 
algorithm (EFF), linear generation mechanism of candidate solution is incorporated with fruit fly algorithm. 
Then status of material algorithm (SMA) is projected to solve the problem in which to establish the evolution 
procedure direction vector operator assigns a direction to every molecule consecutively to guide the particle 
progression. Proposed enhanced fruit fly optimization algorithm (EFF), status of material algorithm (SMA) 
has been tested in standard IEEE 30 bus test system and simulation results show the projected algorithms 
reduced the real power loss considerably. 
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